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ABSTRACT 

Current theories of galaxy formation predict that spiral galaxies are embedded in a 
reservoir of hot gas. This gas is able to cool onto the galaxy replenishing cold gas 
that is consumed by star formation. Estimates of the X-ray luminosity emitted in the 
cooling region suggest a bolometric luminosity of order 10 x 10 41 ergs s _1 in massive 
systems. We have used ROSAT PSPC data to search for extended X-ray emission 
from the haloes of three nearby, massive, late- type galaxies: NGC 2841, NGC 4594 
and NGC 5529. We infer 95% upper limits on the bolometric X-ray luminosities of 
the haloes of NGC 2841, NGC 4594 and NGC 5529 of 0.4, 1.2 and 3.8 xlO 41 ergs s" 1 
respectively. Thus, the true luminosity lies well below the straightforward theoretical 
prediction. We discuss this discrepancy and suggest a number of ways in which the 
theoretical model might be brought into agreement with the observational results. A 
possible solution is that the gravitational potentials of these galaxies' dark matter 
haloes are weaker than assumed in the current model. Alternatively, the present day 
accretion may be substantially less than is required on average to build the disk over 
the Hubble time. Our results are, however, based on only three galaxies, none of which 
is ideal for this kind of study. A larger data set is required to explore this important 
problem further. 
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1 INTRODUCTION 

In models of galaxy formation by hierarchical clustering, the 
disks of spiral galaxies form by late accretion of gas which 
cools from an extended reservoir around the galaxy. In these 
models disks are still expected to be growing at the present 
day. Such models have been very successful in accounting 
for properties of galaxies at optical and infrared wavelengths 
(see, for example, Kauffmann et al. 1993, Cole et al. 1994, 
Baugh, Cole & Frenk 1996, Kauffmann et al. 1999). The pro- 
posed disk formation mechanism should produce a signature 
at X-ray wavelengths which, in principle, is detectable with 
the present generation of satellites. 

At the virial temperature of galaxy halos, the dominant 
cooling mechanism is X-ray Bremsstrahlung. If the cooling 
rate is significant, this flux may be visible as a component 
of diffuse X-ray emission extending well beyond the galaxy's 
optical radius. In this paper, we set out to test this disk 
formation model by looking for X-ray emission from the hy- 



pothesized reservoir of cooling gas around three large, late- 
type galaxies, NGC 2841, NGC 4594 and NGC 5529. 

If we assume that the mass of the present disk has been 
built up by a constant rate of accretion over the age of the 
universe (to), a simple estimate of the X-ray luminosity can 
be made from the binding energy of the disk. 
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where Mdisk is the present mass of the disk, v csc = 
14 [2 ln(r v ir/?"disk) +2] 1 / 2 is the velocity required by mate- 
rial at the edge of the disk to escape the halo (assuming 
a flat rotation curve) and V c is the measured circular ve- 
locity of the disk. Adopting parameters suitable for a large 
spiral galaxy such as NGC 2841 (ie., M disk = 5 x 1O 1O M , 
and V c = 317km/s), this suggests that the halo should emit 
a bolometric luminosity of at least ~ 8 x 10 41 ergs s _1 f^. 
If supernovae and other feedback processes are effective in 
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heating halo gas then even more energy must be radiated 
to assemble the observed disk. Note that eqn. (jl|) predicts a 
minimum Lx oc V c 5 , and that the temperature of the cooling 
gas should scale as the halo virial temperature, Tx oc V^. 
Thus, it is important to search for emission in the most mas- 
sive systems available. 

The picture of hierarchical galaxy formation developed 
by White & Frenk (1991, WF91; see also Kauffmann et al. 
1993, Cole et al., 1994) predicts halo X-ray luminosities more 
directly by computing the cooling rate of g func- 
tion of time. At the present epoch, such models (which are 
tuned to reproduce the observed optical properties of galax- 
ies) predict an X-ray luminosity of about 10 x 10 41 ergs s _1 
for galaxies like NGC 2841. 

To date, studies of X-ray emission from nearby galaxies 
(e.g. Read, Ponman & Strickland, 1997) have concentrated 
on X-ray emission on scales comparable to the optical im- 
age of the galaxy. While extended emission has been de- 
tected in some groups of galaxies (e.g. Mulchaey et al. 1996, 
Trinchieri, Kim & Fabbiano 1994, Trinchieri, Fabbiano & 
Kim 1997), this emission is thought to be associated with 
the group potential. Our aim here is to test for the equiv- 
alent emission in isolated spiral galaxies, which we would 
identify with a cooling flow supplying the galaxy disk with 
gas. A few studies have reported diffuse emission from iso- 
lated elliptical galaxies. However, this has been interpreted 
either as resulting from the expulsion of gas from the galaxy 
(e.g. Read & Ponman 1998, Mathews & Brighenti 1998) or 
as a relic of a collapsed group (Ponman et al. 1994). Fur- 
thermore, since in hierarchical models of galaxy formation 
ellipticals are not necessarily expected to be accreting gas at 
the present day, the existence or otherwise of a diffuse X-ray 
component associated with them does not provide a strong 
test of such theories. For this reason, we have restricted our 
study to isolated spiral galaxies. 

Previous X-ray observations of spirals (see, for example, 
Fabbiano & Juda 1997, Burstein et al. 1997, Pellegrini & 
Ciotti 1998, Brighenti & Matthews 1999) do not place strong 
constraints on the models we aim to test. These studies have 
focussed on measuring X-ray fluxes in regions comparable 
to the optical radius of the galaxy, where the intensity is 
strongest. X-ray emission from galactic stellar sources and 
uncertainties in the spatial X-ray surface brightness profile 
make it difficult to determine what fraction, if any, of the 
observed X-ray flux could be due to a low surface brightness 
cooling flow. 

Since the cooling gas must, at some point, settle into 
the galaxy's disk, some emission is expected from this re- 
gion as the gas cools from the virial temperature of the halo 
to much lower temperatures. For example, in the case of 
NGC 4594 models of the type considered here predict a lu- 
minosity comparable to that measured by Fabbiano & Juda 
(1997), 8.9 ± 1.4 x 10 40 ergs/s in the 0.1-2.4keV band within 
the optical radius of the galaxy. However, likely contamina- 
tion by stellar sources complicates the interpretation of this 
comparison. 

In order to test the models in a manner free from uncer- 
tainties of the type just described it is essential to look for 



assumed H = 50 km s~ 1 Mpc _1 , C = 1, Aq = 0, fi b = 0.06 
and erg = 0.67 



low surface brightness emission at large distances from the 
optical disk of the galaxy. It is also necessary to study the op- 
tically brightest galaxies since the expected X-ray luminosity 
depends so strongly on the circular velocity of the galaxy's 
halo. For this reason, we have chosen for this study nearby 
galaxies with particularly high circular velocities. Most pre- 
vious work has studied relatively low mass galaxies from 
which only relatively weak emission is predicted. 

If diffuse X-ray emission from gas in the haloes of iso- 
lated spiral galaxies proves undetectable, this need not imply 
that the mass of halo gas is small since the gas may simply 
be too diffuse to emit X-rays efficiently. However, it would 
then be unable to cool and flow to the centre of the halo. We 
are thus testing whether the disks of spiral galaxies are be- 
ing built up by the accretion of halo gas at the present day. 
The alternative is that spiral disks were assembled at higher 
redshift and evolved as essentially closed systems thereafter. 

The structure of our paper is as follows. In §2, we search 
for large-scale diffuse emission from massive spiral galaxies. 
We are not able to make any convincing detection; the upper 
limits we establish are almost an order of magnitude below 
the predictions of the WF91 model. In §3 we discuss the 
uncertainties in the models and describe how they might be 
modified to agree with the data. The prospects for future 
investigations are outlined in §4. 



2 OBSERVATIONS 
2.1 Target Selection 

The primary aim of this paper is to search for diffuse X-ray 
emission from nearby isolated spiral galaxies. The galaxy 
must be nearby for the signal to be detectable. Further- 
more, the galaxy must have a high circular velocity in order 
that the energy released by the cooling gas be large. A high 
circular velocity also favours detectability because it implies 
a high characteristic temperature for the emission and so 
less absorption by neutral hydrogen in the Milky Way. In 
addition the efficiency of supernovae feedback is expected 
to decline with increasing circular velocity. 

We created a suitable target list by searching for all 
spiral galaxies with recession velocities below 3000 km s~ 
and with HI 20% line widths greater than 580 km s" 1 . 
Of the sixteen galaxies satisfying these criteria three have 
been observed by the ROSAT PSPC: NGC 2841, NGC 4594 
(the Sombrero galaxy) and NGC 5529. The PSPC data for 
these fields were obtained from the Leicester Database and 
Archive Service (LEDAS). 

• NGC 2841. The Sb galaxy NGC 2841 is the nearest of 
our target galaxies. It has a flat rotation curve extending to 
the limit of HI observations, and is an archetype of galaxies 
with massive dark matter haloes (Kent, 1987). Its star for- 
mation properties have been studied extensively by Young 
et al. (1996) using both Ha and far infra-red indicators. 

• NGC 4594. Also known as M104 or the "Sombrero 
Galaxy," this is a highly luminous Sa seen almost edge on. 
Since this galaxy is bulge dominated, diffuse X-ray emission 
similar to that seen around E or SO galaxies (Forman, Jones 
& Tucker 1985) might be expected. However, X-ray obser- 
vations by the Einstein satellite revealed spectral properties 
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Table 1. Observed properties of the galaxies. Shown for each galaxy are its heliocentric radial velocity, distance, W20 (the width of the 
21cm line at 20% of its peak intensity), TVjj (the column density of neutral hydrogen in the Milky Way in the direction of the galaxy), 
''optical (the 25 B-magnitudes per square arcsecond isophotal radius of each galaxy measured along the major axis), Mb (the absolute 
B-band magnitude) and morphological type. 



Name 


Heliocentric radial 
velocity (km/s)° 


Distance (Mpc) 


w 20 

(km s- 1 ) 11 


N H 

(IO 20 cm- 2 ) e 


^"optical 

(arcmin) d 


M B 


Hubble 

type 


NGC 2841 


638 ±3.0 


13.8 ±5.2 6 


674 ±4 


1.45 


4.1 ±0.04 


-20.6 ±0.2 


Sb 


NGC 4594 


1091 ± 5.1 


24.54 ± 2. 5 C 


753 ±7 


3.67 


4.4 ± 0.04 


-23.0 ±0.2 


Sa 


NGC 5529 


2885 ±5.1 


57.72 ± 2. 5 C 


592 ±5 


1.09 


3.3 ± 0.07 


-21.1 ±0.2 


Sc 



a From the NASA Extragalactic Database (NED). 
6 de Vaucoulers (1979) 

c From heliocentric radial velocity plus peculiar velocity correction assuming H = 50 km s" 1 Mpc" 1 . 
d From the RC3 catalogue (available from NED), corrected for inclination. 
e Dickey & Lockman (1990). 



Table 2. Measured rotation velocity and predicted gas tempera- 
ture for the galaxies in our sample. 



Name 


V c (km s 


— l^a 


T v j r 


(keV) 6 


NGC 2841 


317 


± 


2 


0.287 


± 0.004 


NGC 4594 


358 


± 


4 


0.366 


± 0.009 


NGC 5529 


277 


± 


3 


0.219 


± 0.005 



a Calculated from HI line widths using the expression given in 

Tully & Fouque (1985) 

6 As defined by White & Frenk (1991) 



and core X-ray colours typical of less massive spiral galaxies 
(Kim, Fabbiano & Trinchieri 1992). 

• NGC 5529. This is an Sc galaxy seen almost edge-on. 
Because of its relatively large distance, this galaxy is likely 
to provide weaker constraints than the previous two. 

Distances, line widths, optical radii, and other relevant 
properties of these galaxies are listed in Table [lj The dis- 
tances to two of them (NGC 4594 and NGC 5529) had to 
be determined from their observed heliocentric radial veloc- 
ity; a model of the local peculiar velocity field (Branchini 
et al. 1999) was used to correct the observed velocity to the 
true Hubble velocity using an iterative procedure. Table ^ 
lists the circular velocity of each galaxy as calculated from 
its HI line- width and also the virial temperature of the halo 
estimated, as in WF91. 

2.2 Data Reduction 

Ideally, we would like to compare the X-ray flux from an 
annulus centered on the galaxy (but excluding the optically 
bright region to avoid contamination from stellar emission) 
with similar annuli centred on blank fields. This is not possi- 
ble, however, because different ROSAT fields have different 
particle backgrounds. We have therefore adopted the pro- 
cedure of first subtracting the detected events from a large 
radius annulus in the outer part of the PSPC detector as de- 
scribed in Snowden et al. (1994). Background corrected 'on 
source' and 'off source' fields may then be compared to de- 
termine the flux from the diffuse halo of the target galaxies. 



It is important to note that the background subtraction al- 
gorithm of Snowden et al. tends to produce negative flux in 
the off-source images because of the different point-source 
detection thresholds in the centre and outer parts of the 
ROSAT field. Our analysis takes this effect into account. 

The analysis of the data was carried out using the As- 
terix package. The procedure for analysing each image was 
as follows. Firstly, point sources were detected down to a 3cr 
threshold and excluded by cutting out a circle centred on 
each source. The radius of this circle was set equal to the 
radius of the PSF enclosing 95% of the energy at the posi- 
tion of each source. For each image a background annulus of 
inner radius 0.5° and outer radius 0.7° was then used to es- 
timate the vignetting corrected background count-rate, and 
to create two background-subtracted spectra from source- 
centred annuli with inner radius of 5 arcminutes and outer 
radii of 9 and 18 arcminutes. Choosing a fixed angular radius 
ensures that the blank field analysis is the same for each im- 
age. The outer radius was chosen to include the entire inner 
ring of the PSPC. This is comparable to the predicted X-ray 
emitting region of the closer of our target galaxies (see §3). 
The 9' ring is matched to the X-ray emitting region of the 
more distant galaxy NGC 5529. The innermost aperture was 
excluded from the analysis in order to avoid contamination 
by the diffuse and unresolved stellar X-ray emission of the 
galaxy itself (cf. Read, et al. 1997, for example). The X-ray 
image of each galaxy is compared with the optical image in 
Figure |l[ 

The resulting spectra were analysed with the XSpec 
package to determine the source flux using the mekal X- 
ray emission and wabs absorption models. For each halo we 
adopted a virial temperature of 0.2keV. This is slightly be- 
low the virial temperatures estimated in Table ^, in order 
to obtain conservative estimates of the bolometric luminos- 
ity. The hot halo gas is assumed to have a metallicity of 
0.3Zq (where Zq is the Solar metallicity). Neutral absorp- 
tion column densities varied between the fields and were 
taken from Dickey & Lockman (1990). The only remaining 
free parameter in these models is the overall normalization. 
This procedure appropriately weights the contributions of 
the different energy channels and corrects for the different 
neutral hydrogen cross-sections of the fields. 

The fitted spectrum is then used to calculate the ab- 
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Figure 1. Source images for each galaxy in the sample. The left hand column shows the raw image @w(D(M>@ RglS, rMMFfc^BiO(M)^100iin@gB 
after background subtraction, point source removal and selection of the region from 5 to 18 arcminutes. An optical image of each galaxy 
is superimposed on the X-ray image to give an appreciation of the size of the halo region. 
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Table 3. X-ray fluxes and bolometric luminosities for the haloes of three spiral galaxies. The fluxes and luminosities are determined for 
two annuli, with inner radius of 5 arcminutes and outer radii of 9 (ring 1) and 18 (ring 2) arcminutes. The first column gives the name of 
the galaxy. The flux corrected for blank field offsets is given in the following two columns, and the final two columns list the bolometric 
luminosity of each halo corrected for blank field offsets. 





Bg. Corr. 


Bg. Corr. 


Bolometric 


Bolometric 




Flux (ring 1) 


Flux (ring 2) 


L x (ring 1) 


L x (ring2) 


Galaxy 


(10~ 14 ergs/cm 2 /s) 


(10~ 14 ergs/cm 2 /s) 


(10 41 ergs/s) 


(10 41 ergs/s) 


NGC 2841 


5.62 ± 10.33 


22.90 ± 19.87 


0.03 ±0.06 


0.16 ±0.13 


NGC 4594 


17.96 ± 11.33 


14.21 ± 22.40 


0.44 ±0.28 


0.34 ±0.54 


NGC 5529 


5.91 ± 10.03 


-2.28 ± 19.22 


0.78 ±1.34 


-0.29 ±2.48 



Table 4. The total counts per second in six blank fields. Back- 
ground and point sources were removed in each case and an an- 
nulus from 5 to 9 arcminutes (ring 1) and from 5 to 18 arcminutes 
(ring 2) was selected. 



Field 


Ringl 


Ring 2 


rp201558 


-0.017 ±0.007 


-0.063 ±0.015 


rp700255 


-0.004 ± 0.006 


-0.040 ± 0.013 


rp700319 


-0.021 ±0.005 


-0.024 ±0.011 


rp700376 


-0.007 ±0.005 


-0.003 ± 0.011 


rp700387 


+0.007 ±0.004 


-0.013 ±0.008 


rp700389 


-0.005 ± 0.006 


-0.023 ±0.013 



sorption corrected flux from each halo. These fluxes, shown 
in Table |§| are within the range 0.2 to 2.0 keV in order to 
minimise dependence on the assumed temperature. These 
are converted to pseudo-bolometric luminosities (0.001 to 
10.0 keV) using the conservative temperature of 0.2 keV 
and the distances given in Table [j]. 

For comparison, six off-source fields (Table §) with point 
source targets were selected at random from the ROSAT In- 
ternational X-ray/Optical Survey (RIXOS) (Mittaz et al. 
1998). As we have previously noted, the analysis proce- 
dure tends to produce a small negative count-rate when ap- 
plied to the off-source fields. This is due (at least in part) 
to the different flux limits at which point sources are re- 
moved from the target and background annuli. The mean 
luminosity deficit in the off-source fields was added to the 
measured luminosities of the target galaxies to correct for 
this effect. Alternatively one could add the off-source field 
count rates to the on-source fields and only then fit a model 
spectrum to the resulting counts to determine a corrected 
luminosity. We investigated performing the analysis in both 
ways and found that each procedure gives the same result 
to within 10%. Since the off-source count rates vary from 
field to field it is possible that the true correction required 
for each target field is somewhat different from the mean 
correction from the off-source fields. We therefore measure 
the variance in the off-source field correction and add this 
quantity in quadrature to the errors on the fluxes measured 
in our target fields. The variations from field to field are our 
dominant source of uncertainty. 

The luminosities of the target galaxies are presented 
in Table [| For each galaxy we expect a luminosity around 
10 x 10 41 ergs s _1 . However none of the galaxies is detected 
at a 3cr confidence level in either of the annuli considered. In- 



stead, our data sets 95% upper confidence limits of 0.37, 1.2 
and 3.8 x 10 41 ergs s -1 in the larger aperture for NGC 2841, 
4594 and 5529 respectively. For NGC 5529, the relatively 
weak upper limit is due to the galaxy's larger distance. For 
the two more nearby galaxies, the upper limit lies an order 
of magnitude below the initial theoretical prediction. 



2.3 Uncertainties in the Observations 

In we consider possible modifications to the model of 
WF91 that might alleviate the discrepancy with our data. 
First, we take a closer look at the assumptions we made 
in interpreting the observational data. There are a num- 
ber of uncertain quantities that could affect our conclusions. 
Throughout the analysis we have consistently assumed con- 
servative values for these quantities (i.e. those that would 
minimize the conflict between model and data.) It is never- 
theless instructive to consider how these choices might affect 
our results. 

The dominant source of uncertainty in our analysis is 
the process of background subtraction. Although an accu- 
rate estimation of the background from a mosaic of sur- 
rounding pointings would be preferable, our off-source fields 
indicate the magnitude of the uncertainty. Even using the 
'worst case' field (rp201558), the predicted and observed 
count-rates still disagree by a factor of 2-3. 

Another important source of uncertainty is the extrap- 
olation from ROSAT count-rate to bolometric flux. At the 
low temperatures appropriate to galactic halo gas, neutral 
hydrogen absorption in our own Galaxy has a strong effect. If 
our adopted values for the neutral hydrogen column density 
were inaccurate, this could have a substantial effect on our 
results. However these galaxies are at high Galactic latitude 
and we find that increasing the assumed column density by 
a factor of 3 over the values given by Dickey & Lockman 
(1990) only changes the inferred bolometric luminosity by 
40%. 

The effects of background correction and absorption 
alone cannot reconcile the observed and predicted luminosi- 
ties. Altering the assumed temperature of the halo plasma, 
however, is a more promising alternative. This sensitivity 
arises because, at the relevant temperatures, most of the 
emission lies outside the 0.2-2.0 keV band of the observa- 
tions, forcing us to extrapolate the model fit to lower ener- 
gies. For example, at a temperature of 0.2 keV, 90% of the 
bolometric luminosity is emitted in the ROSAT 0.2-2.0 keV 
band, while at 0.1 keV, only 75% of the emission lies in this 
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band. At even lower temperatures, the effect becomes even 
more pronounced. Applying a lower temperature model to 
the observed spectra, we find that a decrease in temperature 
by a factor of 4.0 (for NGC 2841) to 4.5 (for NGC 4594) is 
required to bring the observations into line with the pre- 
dicted luminosity. It is not clear, however, how these lower 
temperatures can be accommodated within the theoretical 
model. We now turn to a detailed examination of this model. 

3 DISCUSSION 

The model presented by WF91 predicts that the X-ray lumi- 
nosity due to gas cooling within a dark matter halo should 
be: 

Lx = M cool (r cool )V c 2 In ( -^SL J , (2) 

V ^"optical I 

where M coo i(r coo i) is the rate at which mass is cooling in the 
dark matter halo, r coo i is the cooling radius (which grows 
with time), and r op ti C ai is the optical radius of the galaxy, 
which we take to be the extent of the 25 B-magnitudes per 
square arcsecond isophote along the major axis. This model 
assumes that the halo has an isothermal potential with circu- 
lar velocity V c . As shown in the Appendix, the model can be 
readily extended to the case where the gravitational poten- 
tial follows the functional form favoured by Navarro, Frenk 
& White (1996, NFW). 

We have used the formulae and parameter choices in 
the Appendix to predict the bolometric X-ray luminosity 
for each of the three galaxies (using the circular velocities of 
Table ^), as well as the region in which cooling is expected 
to occur. Table ^| lists the cooling radii and luminosities of 
the haloes predicted by WF91's standard model and by a 
model with an NFW potential. In column 5, we list the lu- 
minosity that is expected within the 5-18' annulus used in 
the observational measurements. This geometric correction 
is made under the assumption of an isothermal flow in quasi- 
hydrostatic equilibrium as described in the Appendix. The 
figure given includes a correction for the flux from the target 
galaxy that is expected in the outer background annulus by 
extrapolating beyond the cooling radius with an r~ 2 den- 
sity profile. The corrected flux is about 50% and 70% of the 
total for NGC 2841 and NGC 4594 respectively. The larger 
correction for NGC 2841 arises because the emission from 
the target almost fills the PSPC field of view. The geometric 
correction for NGC 5529 is greater because a large propor- 
tion of the luminosity is projected onto the optical image of 
the galaxy. In the case of this galaxy, the large observational 
uncertainties mean that the emission model is only weakly 
constrained. 

For both NGC 2841 and NGC 4594 the bolometric X- 
ray luminosities predicted by the model are strongly incon- 
sistent with the observations. The WF91 model predicts 
cooling rates of 7.4 and 9.0M Q yr _1 for NGC 2841 and 4594 
respectively. Since the mass cooling rate is proportional to 
the X-ray luminosity in these theories, the simplest inter- 
pretation is that the cooling rates in these galaxies are also 
much less than those predicted by WF91. If we assume that 
the geometric factors would be unchanged by reducing the 
cooling rate, we arrive at upper limits to the present-day 
accretion rate of 0.5 and 1 AAM Qyr~ 1 . 



Further insight can be gained by considering the total 
mass of gas within the cooling radius and the mass of the 
galaxy disk. In the case of NGC 2841 the disk mass can 
be determined by fitting the galaxy rotation curve. Kent 
(1987) derives a disk mass of 9.3 x 10 10 M© allowing a free-fit 
of disk and isothermal halo parameters. Adopting an NFW 
model for the halo, Navarro (1998, N98) derives a best-fit 
disc mass of 5 x 10 10 Mq; a significantly smaller disk mass 
would be inconsistent with the blue luminosity of the disk 
(Young et al., 1996). The WF91 model prediction for the 
mass of gas within r coo i, which is assumed to have formed 
the disk, is 9.8 x 1O 1O M0. Thus the model allows enough 
gas to cool to make the observed disk. Consequently X-ray 
observations imply that the present-day accretion rate does 
not significantly contribute to the total disk mass. Clearly, 
this is an important result, since (if these galaxies are rep- 
resentative) it implies that the masses of spiral disks have 
remained almost constant over recent look-back times. Fi- 
nally, we can compare the gas accretion rate with the star 
formation rate in NGC 2841 of O.8M yr _1 inferred from Ha 
emission (Young et al. 1996). Thus, although our data ex- 
clude an infall rate sufficient to build the disk over a Hubble 
time, they may allow enough infall to replenish the gas being 
used up by star formation. 

Before reaching wide-ranging conclusions, however, we 
must ask whether it is possible to modify the emission model 
to bring it more into line with the observational data. Below 
we consider possible modifications to the halo model that 
might reduce the discrepancy with the observational results. 

Feedback. We have chosen galaxies with particularly 
deep potential wells in order to minimise the uncertainties 
arising from supernovae feedback effects on our predictions. 
According to the models of Kauffmann et al. (1993) and 
Cole et al. (1994) feedback in ~ 300 km/s halos reduces 
the predicted star formation rate by only ~ 5%. Any gas 
ejected from the galaxy by feedback should begin cooling 
once more, and may actually increase the X-ray luminosity 
above our theoretical predictions. Feedback in these models 
cannot, therefore, reconcile theory and observations. How- 
ever, in other models of galaxy formation, such as that of 
Nulsen & Fabian (1995, 1997) and Wu, Nulsen & Fabian 
(1999), in which feedback in earlier generations of dark mat- 
ter halos alters the gas density profile (thereby reducing the 
cooling rate), the expected X-ray luminosities of spiral ha- 
los are likely to be substantially less than in the models we 
have tested. Although these models do reduce the luminosity 
of X-ray halos to levels consistent with our data, the cool- 
ing rate at late times appears much too small to allow the 
growth of disks as large and as massive as those observed. 
Further work is needed to investigate this point further. 

Cosmological Parameters. The predicted X-ray lumi- 
nosities have a relatively weak dependence on the values 
of the cosmological parameters (here we adopted Ho = 
50 km s _1 Mpc _1 , f2 = 1, A = 0, Q, h = 0.06 and erg = 
0.67). The strongest dependence is on the baryon fraction 
since the cooling rate is proportional to the square of the 
density of the X-ray emitting plasma. Using Sib = 0.04 (a 
50% reduction) reduces the luminosities by a factor ~ 2, 
but at the expense of lowering the mass cooling rate also by 
a factor ~ 2. Reasonable changes in cosmological parame- 
ters are unable to lower the predicted X-ray luminosity by 
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more than a factor ~ 3. In particular, lowering Qo tends 
to increase the cooling rate (and consequently the X-ray lu- 
minosity) because the baryon fraction of the halo is then 
increased. 

X-ray Emission Temperature. An assumption of the 
model is that the X-ray emission occurs at a characteristic 
temperature close to that of the galactic halo. If the plasma 
had a multiphase structure, then the emission weighted tem- 
perature could be lower. However, a 25% decrease in tem- 
perature reduces the expected flux by only 20%. Much larger 
changes in temperature are required in order to account for 
the discrepancy. 

Spatial Distribution of X-ray Emission. Because the 
cooling radius estimated for NGC 2841 is so large, the cor- 
rection for the spatial distribution of the X-ray emission is 
important. In our standard model, approximately 60% of 
the total luminosity is emitted within the 5-18 arcminute 
annulus which we observe. If, for example, we considered a 
very different spatial distribution for the luminosity (whilst 
keeping the total luminosity the same) , such as one in which 
the X-ray emissivity remains constant within the cooling ra- 
dius, then the fraction of the total luminosity emitted within 
the annulus would drop to around 30%. The luminosity (cor- 
rected for geometric effects) listed in Table ^ would therefore 
be reduced by a factor of 2. These corrections are much less 
important for NGC 4594 where the cooling radius is better 
matched to the PSPC field of view. 

Halo Circular Velocity. The luminosity of the flow is 
based on the premise that the high circular velocities of the 
disks of these galaxies reflect the circular velocities of the 
halo that confines the hot plasma. Recently, several authors 
(Zaritsky et al. 1997, N98) have suggested that high circular 
velocity galaxies could have formed in significantly cooler 
haloes. N98 has analyzed the rotation curve of NGC 2841, 
and found a best fitting model with V200 = 167j^3 4 km 
s _1 (where V200 is the circular velocity at the virial radius). 
Using the corresponding potential in the calculation of the 
X-ray luminosity we predict a substantially reduced X-ray 
emission of 0.9 x 10 41 ergs s _1 (including geometric factors), 
in much better agreement with the observational results. 
The ratio of the predicted value to the observational upper 
limit is 0.62. Unfortunately such a model underpredicts the 
disk mass inferred by N98 by a factor of 4. Increasing the 
gas fraction to remedy this problem restores the predicted 
X-ray luminosity to an unacceptably high value. 

Metal abundance. We have assumed that the hot gas in 
the haloes has a metal abundance similar to that measured 
in galaxy clusters, namely Z = O.3Z0. Using Z — 0.5Zq 
increases the predicted X-ray luminosity by approximately 
30%, since gas with a higher metal content can cool more 
efficiently. Reducing the abundance to O.IZq would, for the 
same reason, reduce the predicted luminosities (and gas ac- 
cretion rates), by a factor of w 3, but would also reduce disk 
masses by the same factor. 

Remaining gas fraction. The predicted luminosity also 
depends upon the value of / g , the fraction of the baryonic 
material in the halo that is in gaseous form at the time when 
the halo forms. This parameter is unlikely to be much less 
than the value of unity that we have assumed unless a sig- 
nificant fraction of the baryons in galaxies are in some dark 
form. Baugh et al. (1998) find that/ g > 0.7 for halos of circu- 
lar velocity comparable to those of the galaxies considered 



here. If / g were lowered enough so that the predicted X- 
ray luminosity agreed with the observed upper limits, then 
the gas cooling rate would be significantly lower, making it 
impossible to build up a disk of the observed mass by the 
present. Therefore, although we do not know the exact value 
of / g , a lower value cannot reconcile the model with both 
the observed X-ray luminosity and the inferred disk mass. 

Multiphase flow. We have assumed that gas flows into 
the galaxy at the centre of the halo at a fixed tempera- 
ture and only cools once it has settled in the galaxy disk. 
However, it is conceivable that the gas could instead cool 
"in-situ" (i.e. at the location in the dark matter halo where 
it originated) and then fall into the galaxy as cold clouds 
(Thomas, Nulsen & Fabian 1987). If this were indeed the 
case, then the bolometric X-ray luminosity emitted by the 
cooling gas would be 2.5M coo iVq (assuming that the gas 
cools at constant pressure). For typical values of r coo i and 
^optical this is about half the luminosity predicted by eqn. 
(El). Furthermore, this luminosity would be emitted at a 
lower mean temperature than the virial temperature of 
the halo. Using the cooling flow model of Mushotzky & 
Szymkowiak (1988) we estimate that, for gas cooling from 
0.2 keV, 11% of the bolometric luminosity would be emit- 
ted in the 0.2-2.0 keV band after accounting for absorption. 
This is approximately half of the 24% that is emitted in this 
band for gas at a fixed temperature of 0.2 keV. Thus, in 
the multiphase model the predicted flux in the 0.2-2.0 keV 
band is lower by a factor of approximately 4 compared to our 
standard model. These models are thus much less discrepant 
with the observational upper limits. 

It should be noted, however, that even in this scenario the 
infalling gas must still release both its thermal and gravi- 
tational energy somewhere. The discussion above accounts 
only for the thermal component, and this model only suc- 
ceeds if the gravitational component is radiated outside the 
ROSAT energy band. This not readily accomplished: if the 
cold clouds fall ballistically towards the centre, the gravita- 
tional energy would be radiated as post-shock X-rays when 
they collide with the galaxy's gas disk. It does not help if 
the cold clouds are coupled to the hot gas by magnetic fields 
since although they could flow in more slowly by transfer- 
ring their gravitational energy to the hot component, this 
energy will ultimately be radiated in the X-ray band. 
Support for a multiphase model may be provided by the 
high- velocity clouds recently detected by Blitz et al. (1999) 
if these are indeed infalling onto the Milky Way. However, 
Blitz et al. estimate that the present day infall rate of such 
clouds onto the Milky Way is around O.8M yr _1 , and this is 
significantly lower than the rate predicted by the models we 
consider (~ 5-6 Mq yr _1 ). Thus, unless the observed clouds 
represent only a small fraction of the total cloud population, 
the simple multiphase model that we have considered seems 
to conflict with the inferred infall rate of cold gas onto the 
Milky Way. 



4 CONCLUSIONS 

We have searched for large-scale diffuse X-ray emission from 
three nearby, high rotation velocity spiral galaxies. In each 
case, we find that the observed emission is almost an or- 
der of magnitude below that required to build the disk at a 
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Table 5. Predictions from our models, assuming the standard cosmology and halo gas with Z = 0.3Zq. Columns 3 and 4 list the total 
model luminosity and cooling radius in arcmin; column 5 gives the model luminosity after correction for the geometry of the source and 
background annuli; column 6 lists the ratio of the observed luminosity to the model luminosity after correction for the geometry. 



Target 


Halo model 


Model 
L x /10 41 ergs/s 


^cool 

(arcmin) 


Geo. corrected 
Model Lx/10 41 ergs/s 


Ratio of obs. 95% 
upper limit to model 


NGC2841 


Isothermal 


10.4 ±0.10 


35.1 ± 13.7 


5.5 ±0.06 


< 0.07 




NFW 


7.6 ± 0.74 


48.0 ±18.1 


4.5 ±0.44 


< 0.09 


NGC4594 


Isothermal 


10.8 ±0.18 


18.7 ± 1.9 


7.6 ±0.13 


< 0.16 




NFW 


9.5 ±0.63 


23.0 ±2.4 


6.9 ±0.46 


< 0.17 


NGC5529 


Isothermal 


3.4 ±0.03 


9.4 ±0.4 


1.3 ±0.01 


< 2.9 




NFW 


2.4 ±0.10 


12.5 ± 0.5 


0.9 ±0.04 


< 4.2 



constant rate within the age of the universe. The accretion 
rate is similarly discordant with theoretical estimates of the 
cooling rate of gas in an isothermal or NFW-profile halo. 

Although extremely suggestive, our results are not 
definitive. The reasons for this are: 

• Strong constraints are derived only for the two nearer 
galaxies. We may, by chance, have selected two systems that 
have anomalous cooling rates. This is plausible for the bulge 
dominated NGC 4594 galaxy ('the Sombrero'). NGC 2841 
appears to be more typical of spiral galaxies, but its proxim- 
ity makes accurate background subtraction difficult, and our 
results are dependent on the profile assumed for the X-ray 
emission. By contrast, the third galaxy meeting our selection 
criteria is too distant for strong limits on the X-ray flux to 
be set. We note, however, that the X-ray limits established 
for these galaxies are compatible with the emission that is 
thought to occur in the halo of our own galaxy (eg., Moore 
& Davis, 1994, Wolfire, et al. 1995). 

• There are a number of parameters that have been as- 
sumed in order to predict the X-ray flux from the properties 
of the halo. No parameter, or combination of parameters, 
can simultaneously give a low X-ray luminosity and build 
the observed disk. To reconcile the predictions with the ob- 
served luminosities we must assume that either: (a) disks 
were built at high redshift; or, (b) accretion onto disks oc- 
curs without radiating the binding energy of the gas at X-ray 
wavelengths (e.g. in a multiphase flow). 

Our results encourage further work. The current up- 
per limits are set by our ability to subtract the background 
count-rate in the images. This is exacerbated in the ROSAT 
data because the flux level at which individual point sources 
may be subtracted varies across the image on scales com- 
parable to the spiral galaxy haloes. Future work requires 
large-scale imaging with a detector that has a near-uniform 
detection threshold, such as the EPIC camera aboard the 
XMM satellite (Lumb et al. 1996). Such a survey would need 
to target a greater fraction of the nearby galaxies, setting 
firmer limits on the diffuse X-ray emission and eliminating 
variations in galaxy properties as a major source of uncer- 
tainty. If such work continues to fail to make detections, then 
we may be forced to conclude that spiral disks are no longer 
growing at the present epoch. This would agree with tradi- 
tional models of galaxy formation in which galaxies accrete 
their mass at early times and evolve as more or less closed 
systems to the present-day. Clearly, there is much to learn 
from future studies of the X-ray haloes of spiral galaxies. 
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APPENDIX A: A SIMPLE MODEL FOR THE 
X-RAY EMISSION FROM COOLING GAS 

Following WF91, we assume that cooling occurs within a 
static potential, and that the cooling gas remains isothermal 
as it flows to the centre. The assumption of isothermality 
allows considerable simplification in the determination of 
the X-ray luminosity of the flow since it is directly related 
to the drop in potential energy (ie., the internal energy of 
the gas remains constant). Detailed, fully self-consistent flow 
models, such as those of Bertschinger (1989) and Abadi et 
al. (1999), show that these assumptions are reasonable over 
the range of scales over which we wish to determine the flow 
solution. The X-ray luminosity emitted by a unit mass of 
gas as it moves within the potential, <f>, is 



ix = -vr^t- (Ai) 

or 

In models such as those of Kauffmann et al. (1993) and 
Cole et al. (1994) gas is assumed to flow from the cooling 
radius towards the centre of the halo. No mass is assumed 
to fall out of the flow until it reaches the central galaxy. We 
find the total energy emitted by unit mass of gas as it flows 
inwards and then integrate from r coo i to r opt icai to find the 
total energy emitted by the mass of gas, 



r optical 



lx(r) 
Vr(t) 



dr. 



(A2) 



If we now make the assumption that mass flows to the centre 
more quickly than the structure of the flow changes, then the 
structure of the flow at a snapshot in time can be represented 
by the trajectory of a single mass unit. Hence 



Lx = Mcooi(r C ooi) 



dr. 



(A3) 



Below we consider t wo p ossibilities for the structure of 
the dark matter halo. In §A1 



and in 



A2 



we consider an isothermal halo 
we consider the dark matter profile described by 
NFW, in each case assuming that the gas initially traces the 
dark matter. 



Al Cooling in an Isothermal Halo 

This potential was considered by WF91. The dark matter 
density profile is p(r) tx r~ 2 , where r is radial distance in 
the halo. The cooling radius, r CO oi, in the halo is defined as 
the radius at which the cooling time, t coo i, of the gas equals 
the age of the Universe, to- Ideally, we would use the time 
since the halo last doubled in mass (Cole & Lacey, 1993), 
but this information is not available for an individual galaxy. 
Since mass cooling rates decline with time, this leads to a 
conservative comparison with our observational results. 

Assuming a spherically symmetric halo in which the gas 
is initially distributed like the dark matter, the cooling time 
is given by 



i(r) = 



1927rf2 GmpV 



(A4) 



49/ g n b A( A tm p V r c 2 /2k) ' 

(White & Frenk 1991), where tt b is the baryon fraction, A(T) 
is the cooling function of the gas, f s is the fraction of the 
initial baryon density remaining in gaseous form (which we 
normally assume to be 1), G is the gravitational constant, 
m p is the mass of a proton, fi is the mean molecular weight 
of the gas, V c is the circular velocity of the halo and k is 
Boltzmann's constant. 

The evolution of the cooling radius determines the mass 
of gas cooling per unit time. 



M cool {r 

cool J 



1 /gfib K 2 ( ^"cool )f"cool (Vc,max j %) 



2 n 



Gto 



(A5) 



This can be combined with equation (A3) to estimate 
the total X-ray luminosity of the halo. 



A2 Cooling in an NFW Halo 

NFW provide a fitting formula for the density profiles of 
dark matter haloes in N-body simulations of hierarchical 
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clustering. This profile is characterised by two parameters: 
5 C , a characteristic density, and r s , a characteristic radial 
scale in the halo. NFW find that S c and r s are correlated in 
any given cosmology so that the NFW profile provides a uni- 
versal, one-parameter model of the halo. This formula pro- 
vides a good approximation to the density profile of haloes 
of all masses over at least two orders of magnitude in radius 
and its shape is independent of the values of the cosmologi- 
cal parameters. 

Appropriate parameters for the galaxy haloes in any 
given cosmology are determined by requiring that the mea- 
sured circular velocity of the galaxy corresponds to the max- 
imum in the NFW rotation curve. This must be fitted in an 
iterative fashion since the ratio between V CjV ir and V c ,max 
itself depends on c = r2oo / r s ■ 

Using the circular velocity profile of an NFW halo, the 
cooling radius is defined by the relation 

192tt 1 Q Gm^r s 2 a; cool ,2 
*° = — — — y -pr- 7-7= , (1 + CXcool) 

49 fg fib A(r vir ) Xvir 

x [ln(l + cx vir ) - caMr/(l + caw)]" 1 ■ (A6) 
where x = r/r2oo- We then obtain, 

L X = Mcool(r C ool)K 2 ('"200) 

X r COC " X -2 H1 + CX)~CX/(1 + CX) 

L«J Ml + c)-c/(l+c) dX > (A7) 

where 2: CO oi = 7- CO oi/r200 and a; op ticai = r op ticai/»"200- This 
equation must be integrated numerically for any given c. 
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